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Abstract 



It is known that the solutions of pure classical 5D gravity with AdS 5 
asymptotics can describe strongly coupled large N dynamics in a universal 
sector of 4D conformal gauge theories. We show that when the boundary 
Q^ 1 metric is flat we can uniquely specify the solution by the boundary stress 

tensor. We also show that in the Fefferman-Graham coordinates all these so- 
lutions have an integer Taylor series expansion in the radial coordinate (i.e. 
no log terms). Specifying an arbitrary stress tensor can lead to two types 
j> ■ of pathologies, it can either destroy the asymptotic AdS boundary condi- 

tion or it can produce naked singularities. We show that when solutions 
have no net angular momentum, all hydrodynamic stress tensors preserve 
the asymptotic AdS boundary condition, though they may produce naked 
. singularities. We construct solutions corresponding to arbitrary hydrody- 

namic stress tensors in Fefferman-Graham coordinates using a derivative 
expansion. In contrast to Eddington-Finkelstein coordinates here the con- 
straint equations simplify and at each order it is manifestly Lorentz covari- 
ant. The regularity analysis, becomes more elaborate, but we can show that 
^ ' there is a unique hydrodynamic stress tensor which gives us solutions free 

■ of naked singularities. In the process we write down explicit first order solu- 

tions in both Fefferman-Graham and Eddington-Finkelstein coordinates for 
hydrodynamic stress tensors with arbitrary rj/s. Our solutions can describe 
arbitrary (slowly varying) velocity configurations. We point out some field- 
theoretic implications of our general results. 



1 Introduction 

In one of the major developments of late 20-th century physics, it has been shown 
that many strongly coupled conformal 4D gauge theories at large N can be solved 
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by using a classical theory of gravity in ten dimensional spacetime with AdS 5 xX 
asymptotics [1-3]. X is a compact Sasaki-Einstein manifold and is related to the R 
symmetry of the theory if the gauge theory is supersymmetric. In the classical the- 
ory of gravity the dynamics of the metric will be described by Einstein's equation 
sourced by a matter energy-momentum tensor. The matter content of the theory 
of gravity will depend on the presumed dual gauge theory. By the gauge/gravity 
duality any smooth solution of the equations of motion of the classical theory of 
gravity is dual to an on-shell state in the conformal gauge theory and encodes all 
the dynamics of the strongly coupled CFT state in a precise way [2,9]. 

There is, however, always a sector of the theory where the dynamics is uni- 
versal. This is because any two-derivative theory of classical gravity which has 
AdS 5 x X as a solution always admits a consistent truncation to five dimensional 
Einstein's equation with a negative cosmological constant. For instance, we can 
set all scalar fields arising from Kaluza Klein excitations on X and other sources 
to values that minimise the potential and turn off all other matter fields. 

Using AdS/CFT correspondence, now we can define the universal sector of all 
strongly coupled (large N) conformal field theories with gravity duals as follows. 
This sector by definition is the dual of pure 5-dimensional gravity with asymp- 
totic AdS boundary condition. A state in this universal sector will be dual to a 
smooth solution of Einstein's equation with negative cosmological constant. At 
finite temperature also, this correspondence works, but now the solutions of pure 
classical gravity are required to be free of naked singularities Q 

In the first part of the paper we will argue that all solutions of pure classical 
gravity in the universal sector with AdS 5 asymptotics are uniquely determined by 
the boundary stress tensor when the boundary metric is flat. The AdS 5 asymp- 
totics always requires a choice of a boundary conformal structure which means 
that the induced metric on the surface at infinity has a double pole in the radial 
coordinate and its residue can only be fixed upto conformal transformations in the 
boundary coordinates. We say that the boundary metric is flat when we choose 
the boundary conformal structure to be that of flat space. In the gauge/gravity 
dictionary it translates into the dual CFT living in flat space. So our result im- 

'This universal sector is different from what in the context of calculating the tachyon vacuum 
in string field theory is also called the universal sector of 2D CFTs. In the latter case, it is defined 
to be the set of states generated by the action of Virasoro generators on the vacuum [8]. However 
these states cannot be uniquely specified just by the vev of stress tensor alone whereas all solutions 
of pure gravity can be uniquely specified by the boundary stress tensor. So even for 2D CFTs our 
universal sector (which can be defined to be the dual of pure 3D gravity with negative cosmological 
constant) is different from the other definition. 
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plies that in the universal sector the strongly coupled dynamics of the CFT state at 
large N is specified once the conservation of the expectation value of the traceless 
stress tensor is satisfied. From the field-theoretic point of view, this is a surprising 
simplification of the dynamics. 

To establish our claim we will use a theorem due to Fefferman and Graham [4], 
which states that for any solution of Einstein's equations with AdS asymptotics 
we can always use a certain coordinate system within a finite distance from the 
boundary. Skenderis and others [5, 9] have shown that this Fefferman-Graham 
coordinate system also captures the physics of the CFT nicely, in particular, one 
can read off the expectation values of various operators in the dual CFT state and 
also the Weyl anomaly directly from the metric in this coordinate system. We 
will use some characteristics of the CFT to argue that when the boundary metric 
is flat the metric in Fefferman-Graham coordinates should have a simple integer 
Taylor series expansion in the radial coordinate. In fact our argument remains 
valid whenever the Weyl anomaly of the dual CFT vanishes. The result has been 
proved in generality for even dimensional AdS by Fefferman and Graham for any 
choice of boundary metric. Since the Weyl anomaly for any CFT in odd number 
of dimensions vanish, this is a special case of our result. We will use our power 
series ansatz for the metric in Fefferman-Graham coordinates to show that the 
boundary stress tensor expectation value uniquely fixes all the coefficients in the 
power series thus specifying the solution uniquely. Given the CFT argument for 
the consistency of the power series ansatz we will be able to establish that the 
metric is uniquely determined locally by the stress tensor. 

It is clear, however, that any arbitrary traceless and conserved stress tensor 
will not correspond to a CFT state. For AdS 5 asymptotics we can say something 
more about gravity solutions with such boundary stress tensors. Even in these 
cases, we will prove rigorously that the power series solution with no log terms 
in the radial coordinate exists when the boundary metric is fiat. However in such 
gravity solutions either of two distinct pathologies can occur. For stress tensors 
with pathology of the first kind the reverse question of finding the corresponding 
gravity solution will be ill-posed. For such stress tensors, the formal power series 
solution of the metric in Fefferman-Graham coordinates will exist but this power 
series will have zero radius of convergence in the radial coordinate. These patho- 
logical stress tensors will be of the "asymptotic boundary condition destroying", 
or, in short, of "abed" type. The other distinct set of pathological stress tensors 
will produce naked singularities in the bulk. 

We will argue that "abed" type of stress tensors can be avoided by doing a 
perturbation around a stationary late-time solution. We will further specialise to 
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solutions with no net asymptotic angular momentum and these solutions at late 
times will always settle down to a static black brane. H Multi blackbrane static 
solutions will not occur if there are no p-form gauge fields as is the case in pure 
gravity. We will set up a perturbation expansion in the Fefferman-Graham coordi- 
nates and show that all hydrodynamic stress tensors preserve the asymptotic AdS 
boundary condition. This result, we will argue, should also have some measure of 
validity for solutions carrying net angular momentum. 

The perturbation expansion will be similar in spirit to that described in [10, 1 1], 
but we will use Fefferman-Graham coordinates instead of Eddington-Finkelstein 
coordinates. A single black brane preserves the SO(3) rotation symmetries and 
the translation symmetries of the full SO(4,2) isometries of AdS 5 . Among 
the isometries which are broken only two can at most commute with each other 
and there is a four parameter family of choice of these two isometries. Since 
they parametrise the mutually commuting set of broken symmetries of the vac- 
uum, which is the static black brane, we will call these "maximally commuting 
Goldstone parameters". We will choose them to be the scale transformation with 
one scaling parameter and an arbitrary boost parametrised by the three spatial 
components of a velocity. We can use them to generate a four parameter fam- 
ily of so-called boosted black brane solutions. This choice is natural because the 
boundary stress tensor of these boosted black brane solutions will be that of a ho- 
mogenous perfect conformal fluid parameterised by its velocity and temperature. 
The velocity of the fluid will be the same velocity which parametrises the boost 
and the temperature will be the parameter of the scale transformation if the un- 
boosted black brane had temperature unity (in units where the radius of AdS is 
set to unity). Now we will make the velocity and temperature arbitrary functions 
of the field theory coordinates (i.e all coordinates except the radial one) and find 
a correction to the metric which is first order in derivatives of the field theory co- 
ordinates. The boundary stress tensor is also corrected as a result and Einstein's 
equation implies it is conserved and traceless. This perturbation being an order 
by order derivative expansion should be thought of as the holographic dual of the 
usual low energy expansion (E/T) in an effective field theory, T being the temper- 

2 This late time equilibriation, is of course expected only if the boundary stretches indefinitely 
in time, i.e. if the solution is free of "abed" type of pathology. One may see this explicitly by 
studying an example, in which the boundary stress tensor is that of two fluids eternally moving 
past each other at different but constant velocities and temperatures without equilibriating. Our 
results imply that a solution with AdS asymptotics will exist even for such a boundary stress 
tensor. One of the authors (AM) is investigating this solution to check if it indeed has "abed" type 
of pathology. 
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ature. This is therefore a derivative expansion. 

The derivative expansion in the Fefferman-Graham has some advantages over 
the same expansion in Eddington-Finkelstein coordinates [10, 11]. In the Feffer- 
man Graham coordinate system we can naturally view Einstein's equation as evo- 
lution of boundary metric in the radial direction. We will call those components of 
Einstein's equation which contain no more than one derivative of the radial coor- 
dinate as constraint equations. The first advantage is that the constraint equations 
become trivial except for the conservation and tracelessness of the boundary stress 
tensor if the dissipative (i.e the non-equilibrium) part of the boundary stress ten- 
sor t(dis)isv is chosen to satisfy u^t^is^v = 0. The latter is called the Landau gauge 
condition and may be imposed without any loss of generality as by suitable redef- 
initions of the four velocity and temperature we can always make the stress tensor 
satisfy this property. The second advantage over the perturbation in Eddington- 
Finkelstein coordinates is that here the whole procedure will be Lorentz-covariant, 
whereas in the Eddington-Finkelstein coordinates we had to decompose all terms 
into tensors, vectors and scalars of SO(3). The third advantage is that we can 
construct the metric for an arbitrary conformal hydrodynamic stress tensor. We 
can also read off the stress tensor from our metric rather easily. Given this simpli- 
fication of the constraints, in particular, one can think of the Fefferman-Graham 
coordinate system as the "Coulomb gauge" in the context of finding out metrics 
corresponding to arbitrary hydrodynamic stress tensors. 

However, as we already know from the results of [10, 1 1], the solution corre- 
sponding to a generic hydrodynamic stress tensor will contain a naked singularity. 
In the Fefferman-Graham coordinates, however, we will find that the solution al- 
ways has a singularity at the location of the unperturbed horizon. To see if the 
singularity is just a coordinate singularity or a real one we will translate our solu- 
tion to Eddington-Finkelstein coordinates, because in the Eddington-Finkelstein 
coordinate system a real singularity will be manifest in terms of an actual blowup 
of the metric. To do this we will solve the equations of transformation exactly to 
each order in the derivative expansion. We will show that whether the singularity 
in the metric in Fefferman-Graham coordinates is real or fake, the translation to 
Eddington-Finkelstein coordinates can be achieved at every order. The metric in 
the Eddington-Finkelstein coordinates will make the singularity manifest and also 
easily reveal for which choice of the coefficients in the stress tensor would the 

3 The Landau gauge is simply a convenient set of definitions of the velocity and temperature 
variables of the fluid and has nothing to do with gauge fixing of Einstein's equations. The physical 
meaning of these deifinitions is that is the local four-velocity of energy transport. 
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solution be free of naked singularities. At every order in the derivative expansion, 
there will be a unique choice of coefficients of the terms in the stress tensor for 
which the solution will be free of naked singularities. 

Though we will establish the general results stated above, we will give explicit 
computations only upto first order in derivatives. In particular we will find the so- 
lution (exact upto first order in derivatives) in Fefferman-Graham coordinates for 
a conformal hydrodynamic stress tensor with arbitrary rj/s. We will be able to find 
the solution for an arbitrary velocity configuration of the boundary fluid. A special 
case of our result will be the solution corresponding to the Bjorken flow found by 
Janik [13, 14]. With our method we will be able to find the solutions for arbitrary 
slowly varying velocity configurations at each order in the derivative expansion. 
It should also be kept in mind that the pathologies pointed out in [13, 14, 17], 
associated with the methods of finding solutions in Fefferman-Graham coordi- 
nates in [13, 14, 17], do not occur in our case because we never take a late time 
scaling limit in which we are zooming closer to the horizon, where in fact the 
metric always develops a coordinate singularity. In fact our method is as good 
and of equal reach as the derivative expansion in Eddington-Finkelstein coordi- 
nates. It has several comparative advantages which have been pointed out earlier, 
the comparative disadvantage being a slightly more elaborate regularity analysis. 
However if we go beyond the hydrodynamic sector to describe multi black brane 
solutions (if they exist), the Fefferman-Graham coordinate system (being tied to 
the AdS asymptotics) can always be employed efficiently, whereas it is not clear 
if the Eddington-Finkelstein coordinates will be equally useful. 

The organisation of the paper is as follows. In section 2, we establish that the 
boundary stress tensor uniquely specifies a solution of pure classical gravity with 
AdS asymptotics when the boundary metric is flat. In section 3, we confirm our 
claims about the metric in Fefferman-Graham coordinates by translating a known 
solution in Eddington-Finkelstein coordinates which is exact upto first order in 
derivatives and free of naked singularities (we will call this solution as the hy- 
drodynamic solution and has been found in [10]). In section 4, we will set up 
and elucidate the derivative expansion in the Fefferman-Graham coordinates and 
establish that all hydrodynamic stress tensors preserve asymptotic AdS boundary 
condition. In section 5, we will do the regularity analysis of our solutions. Finally 
we will end with some discussion on the field-theoretic implications of our results. 
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2 How the boundary stress tensor fixes the solution 



In this section we will restrict our attention mainly to a five dimensional asymptot- 
ically AdS space with flat boundary metric, though we will indicate in the end that 
our results may be sufficiently generalised. We will soon explain what is meant 
by the boundary metric for asymptotically AdS spaces. 

The Einstein-Hilbert action on 5-dim manifold M, with an appropriate coun- 
terterm to have a well defined variational principle with Dirichlet boundary con- 
dition is 

S=— l — [- f d 5 x^G(R + ^)- f d 4 x^2K] (1) 
lbnO N J M / z J dM 

where K is the extrinsic curvature and y is the induced metric on the boundary. We 
are using the convention of [5] in which the cosmological constant A of AdS d+i 
is normalized to be - d ^\ hence for AdS 5 we have A = -j. 
We want to solve Einstein's equation 

1 __ 6 

Rmn ~ —RGmn - -^Gmn (2) 
1 r 

subject to the condition that the solution is asymptotically AdS with a given con- 
formal structure at the boundary. Fefferman and Graham have shown that for such 
solutions we can use a specific coordinate system called the Fefferman- Graham 
coordinate system near the boundary. In this coordinate system, the metric takes 
the following form: 

ds 2 = G MN dx M dx N = —[dp 2 + gu V (p, z)dz M dz v ] (3) 
P 2 

In the expression above the indices (M,N) run over all AdS coordinates and the 
indices (p., v) run over the four field theory coordinates. The boundary metric 
is defined as 

g(0)iiv(z) = lim^ y (z,p) (4) 

p->0 

Let this boundary metric have a conformal structure. Then it can be shown that any 
conformal transformation of the boundary coordinates (z) can be lifted to a bulk 
diffeomorphism of the Fefferman- Graham coordinates which preserves the form 
of the metric © [6,7]. Under this bulk diffeomorphism, the boundary metric 
undergoes the same conformal transformation. The simplest case for instance 
will be a scale transformation, z — » /tz, of the boundary coordinates for which the 
corresponding bulk diffeomorphism will be p — » Ap (note that in the case of the 
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bulk diffeomorphism, the field theory coordinates z do not transform at all so that 
the boundary metric g mv scales like g (0)MV (z) -» A' 2 g mfJV (z)). 

In the Fefferman- Graham coordinate system the various components of Ein- 
stein's equation reads as [5]: @ 

\g" ~ ~ \g'g~ X g' + -Tr(g- X g')g' - Ric(g) - ±-Tr(g- l g')g = (5) 
2 2p 2 4 2p 

V>( g -y)-v v g;v = () 

Tr[g- l g"] - -Tr[g l g'] - ]-Tr[g- l g'g- l g'] = 
P 2 

Here "(')" denotes a derivative with respect to p and is the covariant derivative 
constructed from the metric g^ v . Also in the above equations we have set our units 
such that 1, the radius of AdS is set to unity. 

When the boundary metric is flat, we will argue that we can expand g MV (z,p) 
in a simple integer power Taylor series of p with coefficients which are functions 
of z. Since we have chosen the boundary metric to be flat, the leading term has to 
be T] M y. Our power series ansatz will be: 

g M v(Z,p) = T]/iv + K=2g(2n)(iv(z)p 2 " ( 6 ) 

We have written down only even powers of p in the above expansion because it 
follows from a result due to Fefferman and Graham [4] that the power series © 
should be an even function of p. H The only even term which is absent is g^viz) 
which follows as an easy consequence of the equations of motion ©. 

It is not obvious that this power series ansatz will indeed provide us a solution, 
so we will give an intuitive argument why this works. By AdS/CFT correspon- 
dence any solution of the bulk equations of motion would give us a state in the 
CFT, so the coefficients of the Taylor series expansion in © should be functions 

4 The (minor) difference with the system of equations given in this reference will be that we 
will use the original Fefferman-Graham radial coordinate p, whereas there the radial coordinate is 
chosen to be the squareroot of ours. Also, the reference uses a definition of the Riemann tensor 
such that the scalar curvature of AdS comes out to be positive. 

5 The existence of power series solution has been proved by Fefferman and Graham for all even 
dimensional asymptotic AdS solutions and in case of odd dimensional asymptotic AdS solutions 
they also argued that if the solution is a power series it should be even. The Fefferman Graham 
coordinates are however unique only upto diffeomorphisms which are the lifts of the boundary 
conformal transformations into the bulk. Although, it is not obvious, it can also be shown [4] that 
the evenness of the series © is independent of the choice of any particular Fefferman-Graham 
coordinate system. 
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of the expectation values of the local operators in the dual CFT state. We will 
explicitly see below that all these coefficients are just functions of the expectation 
value of the stress tensor in the CFT state. It is possible to see the effect of space- 
time independent scale transformation on the CFT operators from g MV (z,p). To do 
this we have to lift the scale transformation to a bulk diffeomorphism so that the 
form of the metric © remains the same and the boundary metric also remains flat. 
This lift, as mentioned before, is achieved by p — » Ap. In the most general case 
it has been shown [9] that the form of the ansatz © should be modified by terms 
like p"(log(p)) m with non-negative n and m. To illustrate our argument we will 
consider just two such possible terms: 

g(n)(z)p n + h (n) (z)p n log(p) 

Under the bulk scaling transformation p — » Ap, 

g (n) (z) -> A n - 2 g in) (z) - log(A)A n - 2 h {n) (z) (7) 

We find the above transformation by checking the new coefficient of p" in g^ v after 
the scale transformation. In a CFT any local operator simply scales like a power 
of A, the power being given by the conformal dimension of the operator. A log(A) 
term is present only when the Weyl anomaly doesn't vanish. In flat space the Weyl 
anomaly vanishes and since we have chosen the boundary metric to be flat the log 
term in © should not be present as g( n )jtv is a function of the expectation values of 
local operators. The absence of the log(A) term in a scale transformation applies 
not only to primary operators but also to their descendents. So we can argue that 
terms like p n (log(p)) m should be absent and g MV should be given by a simple power 
series of p. 

However, our argument, of course, breaks down if the boundary stress tensor 
does not correspond to any CFT state. In Appendix A, we have given the general 
proof of the existence of the power series solution for AdS 5 asymptotics, so that 
even for such cases we can state that the solution, is indeed, a power series. In 
fact we will explicitly see, that for all hydrodynamic stress tensors, whether they 
do or do not correspond to CFT states, the solutions are always power series. 

Now we will substitute our ansatz © in the equations of motion © and solve 
them order by order in powers of p. It is known from earlier work of Skenderis 
et.al. [5] that the first term g^ )llv {z) is just the expectation value of the stress tensor. 
Briefly this is how it comes about to be so. Upto this order the first equation 
(the tensor equation) identically vanishes while the second and third equation of 
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motion give: 

Tr(g (4) ) = (8) 

Since the equations of motion by themselves cannot specify g^ we need a data 
from the CFT to specify it subject to the above constraints. Most naturally g (4) is 
the traceless conserved stress tensor of the CFT. However we can also explicitly 
check this. An explicit calculation shows that g (4) is indeed the Balasubramanian- 
Kraus stress tensor [15] which could be defined for any asymptotically AdS space. 
Hence we may write: 

8(4)fjv = tfjy (9) 

With our ansatz © it turns out that all the other coefficients g^n) (n > 2) are 
fixed uniquely by the equations of motion in terms of g( 4) and its derivatives (or in 
other words the stress tensor and its derivatives). We observe that the first and the 
third of the equations of motion © (i.e. the tensor and the scalar equations) are 
sufficient to solve for g( n y All the higher powers of the second of the equations 
of motion © (i.e the vector equation) identically vanishes on imposing the con- 
straints ([8]) i.e. by imposing the tracelessness and the conservation of the stress 
tensor. It is not difficult to argue that this should be the case because it can be 
shown [5] that the second (i.e the vector) equation of motion simply implies the 
conservation of the Brown- York stress tensor (which when regulated becomes the 
Balasubramanian-Kraus stress tensor) for an arbitrary constant p hypersurface. 
Now the conservation of the Brown- York stress tensor at a given hypersurface is 
not independent of the same requirement for another hypersurface, because in the 
ADM-like formulation of the Einstein's equations if we satisfy our constraints at a 
given hypersurface in which our initial conditions are given the evolution (here in 
the radial coordinate p) automatically satisfies the constraints. The conservation 
of the Brown- York stress tensor at the boundary is already forced at leading order 
in p of the vector equation of motion through ([8]). Hence we should expect that the 
vector equation should not impose any new constraints on the stress tensor given 
that the tensor and scalar equations specify all the coefficients uniquely and this is 
exactly what is borne out. In our proof in Appendix A, we show how the tensor, 
vector and scalar equations of motion turn out to be consistent with each other 
when we employ the power series ansatz. 

Below we give the a few of the the coefficients g^y 

1 
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- -^^(d/A^y + dydjfr) 
oU loU 

+ ^oy^/jv + a v f p d a t Pll ) 



g(12*iv = - ^t^tap) + (10) 

Here, as before in ([5]) the boundary indices are raised and lowered by t]^ and □ is 
the Laplacian in flat space. Let us observe and explain certain simple features of 
the results above. The first observation is that every term in the RHS of the above 
equations contain only even number of derivatives. This is so because the terms 
containing derivatives originate only from Ric(g) in the first of the equations of 
©. The second observation is that the terms independent of the derivatives appear 
only for g^„). This is so because if we omit Ric(g) in the first of the equations of 
([5]), then the solution is a power series in p An as the first non-trivial term in the 
series is g (4 ). So for a solution where the stress tensor is uniform (like in the case 
of a static black brane solution), g has an expansion containing only p 4 " terms. 

With our argument that the ansatz © should give us a consistent solution, it 
is obvious that the stress tensor, which appears as g(4) in g uniquely specifies the 
solution because all the higher coefficients are fixed uniquely in terms of with 
no new constraints like ([8]) appearing for g (4 ). This completes the argument that 
when the boundary metric is flat we should have a solution uniquely specified 
locally by the stress tensor alone. This statement readily generalizes to other 
dimensions in the case of a flat boundary metric and most likely also generalizes 
when the boundary metric is not flat. The general validity could be argued for on 
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the basis of the equations of motion © which are second order (specifically in 
derivatives of p). Intuitively the boundary metric and the stress tensor specifies all 
the initial data we need for a unique solution, however a concrete demonstration 
of this would probably require methods beyond what we have employed here. 

Our power series ansatz © should work even if the Einstein-Hilbert action 
with negative cosmological constant receives higher derivative corrections pro- 
vided the boundary stress tensor corresponds to a state in the dual theory. Our 
argument as to why it works is independent of the equation of motion and like- 
wise also independent of say, the value of t'hooft coupling of the dual theory. We 
have just used the fact that a conformal transformation in the boundary should 
have an appropriate lift to a bulk diffeomorphism consistent with the transforma- 
tion of CFT operators. The transformation of the CFT operators under conformal 
transformations, as well, is independent of the value of the coupling. In fact one 
can readily check that exact static black hole solutions of Gauss-Bonnet gravity 
which are asymptotically AdS (given in [19]) have power series expansion when 
written in Fefferman- Graham coordinates. 

The argument we have given above, however, cannot be reversed to argue 
that a solution with asymptotic AdS 5 boundary conditions exists for any arbitrary 
stress tensor. The reason that we can't reverse the argument is that the series © 
for g^y exists only formally. The coefficients g w may not be well behaved at large 
n, for an arbitrary stress tensor. We will give a simple example to show what 
can go wrong. For a specific choice of stress tensor, we may find that g {n)fJV = 
f( n )s^v plus other terms. Here s MV is a specific term in the stress tensor. If, for 
instance, the series 2,„f(n)p n has zero radius of convergence, g^ v will not be a 
meaningful series of p as it will also have zero radius of convergence in p. Such 
boundary stress tensors, for which g^ v has zero radius of convergence in p, could 
be appropriately called, "asymptotic boundary condition destroying" stress tensor 
or in short "abed" stress tensor. We will have more to say about such stress tensors 
in section 4@ 

interestingly, Fefferman and Graham have shown in [4] that for even dimensional asymptotic 
AdS solutions, always has a finite radius of convergence in p. However their argument does 
not readily generalize to the odd dimensional case. 
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3 Mutual translation between Eddington-Finkelstein 
and Fefferman-Graham coordinates 



In the previous section, we have seen that, the Fefferman-Graham coordinate sys- 
tem is good for finding a solution to Einstein's equation with a negative cosmolog- 
ical constant when the corresponding boundary stress tensor is specified. How- 
ever the solutions are usually found in other coordinate systems. For instance, the 
static black brane solution is usually described in the Schwarzchild-like coordi- 
nate system and the hydrodynamic metric of [10] has been found in the Eddington- 
Finkelstein coordinate system. It would be useful to see how we can rewrite these 
solutions in the Fefferman-Graham coordinate system asymptotically. We will 
demonstrate a novel technique towards this end for the boosted black brane and 
the hydrodynamic metrics. In both cases we will see that we can achieve a mu- 
tual translation between Eddington-Finkelstein coordinate system and Fefferman- 
Graham coordinate system by using a power series ansatz similar to © and we 
can solve this ansatz algebraically order by order. We expect this method to work 
for all solutions in which the boundary metric is flat, or more generally when the 
Weyl anomaly vanishes. 

The general procedure is as follows. In the Eddington-Finkelstein coordinates 
(jf, r) the metric takes the form: 

ds 2 = -2u^ l {x)dxf I dr + G^ v {x, r)dx IJ dx v (11) 

Here we are using ingoing Eddington-Finkelstein coordinate system, so that u M is a 
four- velocity (hence u ll u v rf v = -1) such that it is directed forward in time. We will 
express the general structure of coordinate transformation from the Eddington- 
Finkelstein coordinates (xf, r) to Fefferman-Graham coordinates (z M ,p) as below: 

dp = p M (r, x)dx M + q(r, x)dr (12) 
dz M = m M v (r, x)dx v + n\r, x)dr (13) 

We substitute the above in the Fefferman-Graham form of the metric © to get: 

ds 1 = -^[(PfiPv + gtfip, z)m J] il m\)dx* l dx v + lip^q + g^(p, z)m ' n^dx^dr 
+ (q 2 + g^zh^dr 2 ] 
Comparing the above with the Eddington-Finkelstein form of the metric (fTT|) . we 
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get the following set of equations: 



(q(x, r)f + g MV (p, z)rf(x, r)n v (x, r) = (15) 
r)q{x, r) + g a/3 (p, z){m a ^(x, r)nP(x, r) + rrf^x, r)n a {x, r)) 



= -lu^ipix, r)) 2 



p^x, r)p v (x, r) + g a p(p, z)m ff jU (x, r)rrt \{x, r) = G^ v (x, r)(p(x, r)f 



So we have a scalar, a vector and a tensor equation and three unknowns to solve 
for. The unknowns are a scalar p(x, r), a vector z^ix, r) and the tensor g^ v {z,p) 
which appear in the Fefferman-Graham metric ©. It is clear from the definitions 
(PT2"1) of q, etc. that they are just various partial derivates of (p, z), for instance q = 
d r p, etc. We will make the following general ansatz to solve the above equations. 
The ansatz for p and zf will be that they will be an integer power series of the 
inverse of the Eddington-Finkelstein radial coordinate r. 

1 piix) P3(X) 

P=- + ^- + ^^ + ( 16 ) 

r r z r 3 

zl(x) f 2 (x) 

f = x? + - — + -At- + 

r r l 

To solve the equations of transformation (IT3T ), the above should be supplemented 
with the ansatz © for the gfi V (z,p) in the Fefferman Graham metric. The expres- 
sions for the partial derivatives like q, etc. then turn out to be as below: 

1 2p 2 3p 3 
q = d r p = -— . — - .... (17) 

Y' y 

P M = d M p = — + — + 

r r 1 

One thing to be kept in mind is that when we substitute our ansatz (fT6l ) to solve 
the equations of transformation (fT5l) , g^ip, z) should be re-expressed as functions 
of (x,r). Below, we just give the first three terms which appear after it is rewritten 
as functions of (x,r). 

tftvix) (4p 2 t MV + (zi.d)t MV )(x) 
g MV = %y + — — + ~ 5 + .... (18) 
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We now consider a boosted black brane metric in Eddington-Finkelstein coordi- 
nate 

„2 _ o,, j,jij„ ,2-c/i — \,. ,. j,jij„v . „2i 



where 



J* = -lu^dx^dr - r f(br)u IJ u v dx M dx y + r P fIv dx M dx v (19) 



m = i - ^ (20) 

«° = -=L= (21) 

Pi 



(22) 



and the temperature is T = and the three-velocity /?,- are all constants, and 

*V = «//"v + i7/w (23) 

is the projector onto the spatial hypersurface orthogonal to the four velocity vt. 
This metric can be obtained by applying a boost parameterised by the three- 
velocity Pi and a scaling by b to the usual AdS black hole with unit temperature 
where the time coordinate t is itself a Killing vector. In this case actually the ex- 
act transformation from Eddington-Finkelstein to Fefferman-Graham coordinate 
system can be exactly worked out easily and it is given by: 

, - ■ ^ (24) 
yjb 2 r 2 + VfcV - 1 
? = x? + u^bk(br), 

1 y + 1 

&O0 = t(w£( r) - larctaniy) + n) 

4 y — 1 

The solution for g pv in the Fefferman-Graham metric © for the boosted black 
brane is given by: 

p 4 4p 4 
Spvfcp) = (1 + ^tyv + 4M+p 4 M ^v (25) 

The boundary stress tensor could be easily read off by looking at the coefficient 
of p 4 after Taylor expanding the RHS of the above expression. The stress tensor 
turns out to be that of an ideal conformal fluid (like that of a gas of photons) 

%v = g(4)^ v = -^[Au^Uy + 7fc v ] (26) 
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where the temperature is T = -4. The horizon in the Fefferman-Graham coor- 
dinates is at p = y/lb and at the horizon g MV given by (1231 ) is not invertible as 
8tiv(P = V2Z7, z) = IPpy. So clearly the Fefferman-Graham coordinate system has 
a coordinate singularity at the horizon. Also it is easy to check from (|24l) that the 
change of coordinates also becomes singular at the horizon. 

Now we turn to the hydrodynamic metric found in [10] which is a solution to 
Einstein's equation upto first order in the derivative expansion and has a regular 
horizon. Here the "maximally commuting Goldstone parameters" of the boosted 
black brane solution, the velocities /?' and the temperature T are functions of the 
field theory coordinates (x). The in the Eddington-Finkelstein form of the 
metric (fTTT) is: 

1 2 

G^y = r 2 P^ v + (-r 2 + -g^u^Uy + 2r 2 bF{br)a^ v - r((u.d)u ft Uy - -u^Uyid.u)) (27) 

with 2 

F(x) = -(log(- ^— -) - 2arctan(x) + n) (28) 

4 x 4 

In this case we will solve the set of equation (fT51) by putting in our anstaz (fT6l) . We 
solve order by order for each power n in r~ n . At each order we have to solve alge- 
braic equations and remarkably the equations can be consistently solved at each 
order. It is important to throw away all the terms which have two x-derivatives or 
more and solve the series for p and z M given in (fT6l) and the series for g MV given 
in © only upto first derivative order. This is justified because the hydrodynamic 
metric above in Eddington-Finkelstein form is a solution to Einstein's equation 
only upto first order in x-derivatives and hence it can have a Fefferman-Graham 
expansion near the boundary only upto first derivative order. The results of the 
non- vanishing terms in the expansion for p and z M in (fT6l) upto r~ 9 order are given 
below: 

1 1 I3(d.u) 7 

P2 = ^(d.u), P5 = —, P6 = - T ^ r ,p 9 = — (29) 

1 u? 
z\ = uP,z% = -u> J (d.u),zf 1 5 = jg, 

„ _ 9u^(d.u) + 7(a.y ji _ u 11 
6 ~ 60b 4 9 ~ 9b* 

We can easily observe some patterns in the results above. Firstly the terms with- 
out any derivatives only appear as coefficients of r" 4 " -1 . These are precisely the 
terms that appear in the expansion for the case of the boosted black brane as given 
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in (1241) . This is because the original black brane solution in Fefferman- Graham 
coordinates as we know from (T25T ) is a series with "gaps" of four (which means 
only the fourth next term is non-zero). So the solution of (fT5l) should provide a 
series for p and z M in gaps of four as well. Secondly, it also turns out that the terms 
which have first derivative pieces occur for p2,P6> z£> etc. again in gaps of four. 
We obtain the coefficients of the series for g^ given in © which was part of our 
ansatz. The second non-zero term in the series gives us the boudary stress tensor: 

_ _ V + *u M Uy 1 

fpv - g{A)»v (iU) 

where 

<V = P, tt Pfd (aUji) - l -P, v d a u a (31) 

This is stress tensor for a relativistic conformal fluid satisfying Navier-Stokes' 
equation and with rj/s = 1 1 An. The next non vanishing term in the series for g^ v 
is: 

We can check that the expression for is given by the general results of the the 
previous section when we substitute the dissipative stress tensor (13"0"1 ) in (fTOT ) . 

In this section we have worked out the case for a specific "hydrodynamic met- 
ric" given in [10]. This metric has no naked singularities and this corresponds 
to the choice of tj/s = l/4n in the dissipative stress tensor (|3T| ). However we 
will see in section 5 that our ansatz (fT6l) for translation between the Eddington- 
Finkelstein and Fefferman-Graham coordinates will work even when the above is 
not the case, i.e the metric contains naked singularities. In what follows we will re- 
verse the translation. That is, we will work out the Fefferman-Graham form of the 
metric exactly upto first order in derivatives first and then find out the Eddington- 
Finkelstein form of the metric also exactly upto first order in derivatives. We will 
see that the power series ansatz (fT6l) is consistent for any metric corresponding to 
an arbitrary hydrodynamic stress tensor. 



4 The derivative expansion in Fefferman-Graham 
coordinates 

We have already seen that the Fefferman-Graham form of the metric is the ideal 
one to use if we are asking given a boundary stress tensor what the corresponding 
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solution of Einstein's equations of motion should be. The most general hydrody- 
namic stress tensor for a conformal fluid (in the Landau gauge) upto first order in 
derivatives is as below: 



V + 4m p (z)m v (z) 
Ab{zf 



2b(z) 



7 



(33) 



with cr^yiz) given by (|3TT) . b related to the temperature through b = 1/nT and y 
an arbitrary constant. However here, unlike in the case of the specific solution 
(without naked singularities) we considered in the previous section, rj/s = y/4n 
and hence is arbitrary. We now ask what would be the corresponding solution for 
this arbitrary case. 

Before we get into this specific case, we will show that we can get some in- 
sights into the reverse question from some generally known facts and our previous 
results given in section 2. We have seen, briefly, at the end of section 2 that the re- 
verse question is ill posed for an "abed" (asymptotic boundary condition destroy- 
ing) stress tensor, for which the formal power series © for g MV has zero radius of 
convergence in p. One must devise a strategy in which such stress tensors do not 
appear at all. To this end we may always exploit a general property of solutions 
of Einstein's equation that in the long run the solution always becomes stationary. 
For the moment let us further restrict to those solutions which have no (ADM) 
angular momentum or any other (ADM) conserved charges (like the R-charge). 
These will, in the long run, settle down to the known boosted black brane solution 
(PT91 ). Static multi blackbrane like solutions do not appear if we turn off p-form 
gauge fields, so if more than one black brane are present they eventually will col- 
lapse to form a single black brane. A good strategy to recover all solutions will 
be to perturb around the late-time static black brane and build up all solutions in 
a systematic derivative expansion. Since any solution would eventually become 
static (or equilibriate) this strategy should always work at sufficiently late times. 

Since the approach to equilibrium can be naturally described by hydrodynam- 
ics, one can intuitively expect that the late time behaviour of the solutions will 
correspond to a hydrodynamic description in terms of the boundary theory if the 
equilibrium can be described in terms of a perfect fluid. The boundary stress 
tensor of a boosted black brane indeed corresponds to that of a perfect confor- 
mal fluid like that of photons in pure QED. Our expectation is indeed borne out 
by the fact that all solutions in the derivative expansion correspond to a trace- 
less conserved hydrodynamic boundary stress tensor, but with arbitrary number 
of derivatives. We will see that in the derivative expansion at each order the solu- 
tions always have finite radius of convergence away from the boundary, so we can 
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conclude that all hydrodynamic stress tensors are asymptotic boundary condition 
preserving. 

The fact that all hydrodynamic stress tensors preserve the asymptotic AdS 
boundary condition should have a certain measure of validity even for solutions 
with net angular momentum. In fact in [18], it has been shown that a large class 
of rotating black holes in AdS can be described by perfect fluid hydrodynamics. 
However, we do not know how general the result is. The argument in the previous 
paragraph shows that for any solution if the hydrodynamic description holds for 
the stationary solution to which a given solution eventually equilibrates, it should 
hold for sufficiently late times as well. So certainly a large class of solutions even 
in the sector with net angular momentum which can be constructed by perturbing 
around certain stationary solutions will have a hydrodynamic description at least 
at late times. 

To build up a solution corresponding to an arbitrary hydrodynamic stress ten- 
sor, we will work in the Fefferman- Graham coordinate system as we have said 
before and we will construct the solution exactly order by order in the derivative 
expansion. To develop the derivative expansion we follow the same method which 
the authors of [10] followed but now in the Fefferman- Graham coordinate system. 
In fact, based on the results of section 2, we will see that their method simplifies 
in these coordinates. We take the boosted black brane solution with g MV of the 
form of (|25l ), but now the "maximally commuting Goldstone parameters" {uf 1 , b) 
are arbitrary functions of z. We will call this the zeroth order metric g which is 
no more a solution to Einstein's equation, so we need to correct this with g\ which 
will now depend on the first derivatives of the "maximally commuting Goldstone 
parameters" (u M , b). This correction g\ can be found substituting g = g Q + g\ in 
our equations of motion © and retaining only terms which have no more than 
one derivative of z. 

The first of the equations of motion ©, i.e the tensor equation gives us a 
source free linear equation for g\ which is second order in the derivatives of p and 
has no z-derivatives. 

\g[ ~ I— - ~g'i£oYo - Vi + rgigoWoYo 

P , , (34) 

+ \{Tr(gJg\) - Tr{g 'g Y g- l g Q )){ g -± - *°) + l -Tr{g- Q l g )(^ - *1) = 

7 As we have mentioned in a previous footnote, a non-trivial check of this strategy will be to 
construct a solution for a boundary stress tensor for which there is no late time equilibriation and 
see how it is connected to the "abed" type of pathology. 
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At the first order in derivative expansion, the only term which can provide a source 
term is Ric(g) since it has no derivatives of p. However Ric(g) contains at least 
two derivatives of z, so at this order the source vanishes. 

At the first order the second of the equations of motion, which is a vector 
equation gives us the following: 

V«*,7Xfc Vo) - V V o8o» v = (35) 

where V is the covariant derivative constructed from go-The major simplification 
which occurs in the Fefferman-Graham coordinates is the general observation in 
section 2, that this gives us nothing but the conservation of the stress tensor. It 
may be checked that if we choose to solve this vector fluctuation equation order 
by order in powers of p, like we did in section 2, at the leading order we would get 
d^topy = 0, where ?o//v is the perfect fluid stress tensor (l26l) and all the coefficients 
of the higher powers of p will vanish identically once the leading order condi- 
tion is imposed. This simplification will happen at every order in the derivative 
expansion, which means that if t n -\ is the stress tensor upto n-1 th order in the 
derivative expansion, at the n-th order the second equation will simply imply the 
conservation of t n -\. 

At the first order in the derivative expansion the third equation of motion van- 
ishes identically. It is easy to see why this will happen. Again we go back to 
the general observations of section 2. If t MV = t 0MV + t\^ y with t 0fIV given by the 
perfect fluid stress tensor (|26l) and t lMV is the first order correction to the stress 
tensor satisfying the tracelessness and the Landau gauge u^t^y = conditions, 
then the correction to the coefficients of the power series expansion g( n)MV (some 
of which are listed in (flOl) ) is simply proportional to ti MV . The first order deriva- 
tives of to^ doesn't appear because, as we have observed the general expressions 
for g(n) must contain even number of derivatives of It follows that the cor- 
rection to the zeroth order metric, gi, is proportional to t\. It also follows from 
the the tracelessness of t\ and the Landau gauge condition that the third equation 
vanishes identically as all traces appearing in the equation vanish. We will soon 
see that, this simplifying feature also, remarkably generalises to all orders in the 
derivative expansion. 

In the Fefferman-Graham coordinates the first order correction to the metric 
gi is, therefore, proportional to the first order correction to the stress tensor which 
is proportional to cr^ and therefore gi takes the form of y bcr^fip), where y is 
an arbitrary constant. Substituting this in the tensor equation (fJH), we find that 
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f{p) satisfies the following differential equation: 



„ , (12fr 4 - p 4 )(4b* + 3p 4 ) 128p 6 fr 4 _ 

1 1 p(16Z> 8 -p 8 ) 1 (4b 4 + p 4 )(l6b s - p s ) { ) 

We already know that the solution is a power series in p 4 , so we change our vari- 
able p to x = p 4 . The equation now reads 

_ ■ *# + f 

J 7 16£ 8 -* 2 J (4M + x)(16& 8 -jc 2 ) V ; 

The solution of this differential equation which vanishes at the boundary (after 
resubstituting x with p 4 )[§ is: 

(1 + W )h8( T^ ) (38) 



The metric in Fefferman- Graham coordinates upto first order then is: 

ds 2 = 



dp 2 + g MV (p, z)dz M dz v 



P 2 

p 4 4p 4 / p 4 1 - ^7 

gflv (p, z) = (l + — ) v + — - u^Kv + r'^v(i + -7-) (39) 

To read off the stress tensor upto first order, we simply need the p 4 term in the 
Taylor expansion of g^ v . We get: 



H + 4^ v _ 



Comparing with (1331) we get that we must set y = y in the first order metric ([391) 
to get the desired solution corresponding to the boundary stress tensor. 

One very interesting feature of our solution at the first order can be found 
out by putting y = y = 0. This implies that our zeroth order solution itself, 
now with velocities and temperatures satisfying the relativistic Euler equation, 
is an exact solution of Einstein's equations upto first order. Such is never the 
case in Eddington-Finkelstein coordinate system where as we will see we need to 



8 The other solution is fa — 1 + ^ 
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correct the zeroth order solution even for a dissipation-less stress tensor so that 
the solution is exact upto first order. We do not understand any deep reason for 
this feature of our solution. 

Now we can proceed to examine the higher orders in the derivative expansion. 
Though we will postpone explicit solutions beyond the first order for a future 
publication, here we will show that it is trivial to satisfy the vector and scalar 
constraints at each order in perturbation theory. The tensor equation takes the 
following form at each order in perturbation theory: 

Dyg nliv + D 2 (g nMp U p Uy + g nV pU P U^) + D 3 (g npa .Tf (r )j] flv + D^g^rf^U^Uy 

+ D 5 (g npa -u p u°~)r] flv + D 6 (g npa ifii r )u ll Uy = s nftv (z,p) 

where D\, D 2 , etc. are linear differential operators involving derivatives in the 
radial coordinate only and s nfiy (z,p) is the source term which is a (nonlinear) func- 
tion of the corrections to the metric upto n-1 th order in the derivative expansion. 
The left hand side of the above equation is in fact the same as in (1341) with g\ 
replaced by the n-th order correction to the metric g n , but now source terms are 
present on the right hand side. Also the differential operator Di is the same as the 
operator which acts on f in (|36l) at every order in the derivative expansion. We 
dropped the operators D 2 , D 3 , etc. at the first order, i.e. for gi, because as we 
saw the general results of section 2 (equations in (PTOl) for instance) forced it to 
be proportional to be stress tensor and hence be traceless and vanish when con- 
tracted with the four velocity. However, from the second order in the derivative 
expansion onwards, the general results of section 2 do not imply this to be true for 
the correction to the metric and in fact the source terms which appear on the right 
hand side of the equation indeed do not have this property. All the other operators 
except Di, however, involve no more than one derivative in the radial coordinate. 

We have to choose a particular solution to the above equation. We can always 
choose the particular solution to be such that it vanishes at the boundary like p 6 
so that it doesn't contribute to the stress tensor (as the coefficient of its p 4 term 
vanishes). One can explicitly check this, however, more efficiently we can prove 
it as follows. The source term for the n-th order correction clearly is determined by 
various terms of the stress tensor upto n-1 th order, so it follows from the general 
results of section 2 that the particular solution can be chosen to be independent 
of t nflv , which is the n-th order correction to the stress tensor. In that case the p 4 
term should be absent. For instance, based on the results like those in (PTOl) . we can 
write down the Taylor series expansion in the radial coordinate for the particular 
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solution for g 2 as below. 



P 6 



1 



+ P w [-—(t "pt 0av + t °nt 0ail ) 

1 nR 1 



+ dydJofr) (42) 
+ -^t^dadptQ^y - Y^d^dytoa/s + j^T] fIv d a t^ > 'd a tQ/ 3 y 

+ -^id^datopy + dyt^datopfj) - -^dat^dpt^] + 

More generally, the particular solution for g n is uniquely determined once we 
specify that it vanishes at the boundary like -(l/12)p 6 n?„^ 2 - Then it follows that 
it is independent of t n and doesn't contribute to the stress tensor at the n th order. 

Now the particular solution at every order in the derivative expansion should 
by itself satisfy the scalar constraint. Let us see it explicitly for the particular 
solution for g 2 . The particular solution chosen to vanish at the boundary like 
-(l/12)p 6 n?o has an expansion of the above form (42). So by this choice, the 
coefficients of the Taylor expansion (now fixed by the source) will automatically 
agree with the general formulae, like those in (PTOl) . These general formulae are 
automatically consistent with the scalar constraint. The scalar constraint also will 
be a linear differential equation for g n with a source term. The source term again 
is a (nonlinear) function of the corrections to the metric upto n-1 th order in the 
derivative expansion. The particular solution by itself will satisfy this equation. So 
the homogenous solution of the tensor equation for g n must also be a homogenous 
solution of the scalar constraint. 

The homogenous solution of the tensor equation for g n which will be consis- 
tent with the scalar constraint is simply -2b 4 f(p)t n ^ v , with f(p) being given by 
(TJHT ) and t, vv being an arbitrarily chosen correction to the hydrodynamic stress 
tensor involving n derivatives of the field theory coordinates z. However t nfiv must 
be traceless and also satisfy the Landau gauge condition. Let us illustrate again 
by explicitly doing the Taylor series expansion of the homogenous solution to g 2 
which is -2b 4 f(p)t tltlv . The Taylor expansion is as below: 

8 12 

g2, V = + £j + ^ + -..) (43) 

Using the tracelessness and Landau gauge condition for t 2 , one can check from the 
general formulae like those in (flOl) that this is just the part of the metric determined 
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by t 2 at the second order. Hence this should be the only homogenous solution that 
is consistent with the scalar constraint. Similarly at each order one can see that 
the part of the solution for g n which contains t n is proportional to t n and since the 
particular solution by choice contains all other terms, the homogenous solution 
should be always proportional to t n . Then the tensor equation fixes the radial part 
of the homogeneous solution so that it should be -2b 4 f(p)t, vv . 

The vector constraint, at the n-th order in the derivative expansion, as we have 
argued before simply implies the conservation of the stress tensor upto n-1 th 
order. 

To summarize, these are the features of the derivative expansion in the Feffer- 
man Graham coordinates. 

• At every order in the derivative expansion, the tensor equation for g„ is a 
linear differential equation of the form of (I4TT) involving derivatives in the 
radial coordinate only. The operators Di,D 2 , etc are the same at every order, 
while the source term s n is a nonlinear function of the various corrections to 
the metric upto n- 1 th order. 

• The particular solution to the tensor equation for g n can be chosen to vanish 
at the boundary like -(l/12)p 6 n?„_ 2 . With this choice the particular solution 
automatically satisfies the scalar constraint. 

• The homogenous solution to the tensor equation which is consistent with 
the scalar constraint is -2b 4 f(p)t nMV at very order, with f being given by 
(|38l ) and t nfiv being an arbitrary n th order correction to the stress tensor 
which satisfy the tracelessness and the Landau gauge condition conditions. 

• The vector constraint at the n-th order just implies the conservation of n-1 
th order stress tensor. 

• We can keep manifest Lorentz covariance at each order in the derivative 
expansion. 

• We can construct a solution corresponding to an arbitrary stress tensor be- 
cause the homogenous solution of the tensor equation for g n at the n-th 
order is simply proportional to an arbitrarily chosen n-th order correction to 
the stress tensor. At every order in the derivative expansion for any choice 
of the hydrodynamic stress tensor, the solution has finite radius of conver- 
gence away from the boundary, so all hydrodynamic stress tensors preserve 
the asymptotic AdS boundary condition. 
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5 Getting rid of naked singularities 



The comparative advantage of solving Einstein's equation of pure gravity in Fef- 
ferman Graham coordinates in the derivative expansion over doing the same in 
Eddington-Finkelstein coordinate system is that the constraints simplify dramati- 
cally and also we do not need to split the terms into tensors, vectors and scalars 
of SO(3), thus preserving manifest Lorentz covariance. The comparative disad- 
vantage of the Fefferman-Graham coordinate system is that the regularity analysis 
is not straightforward. At the first order in the derivative expansion, the metric 
in Fefferman-Graham coordinates (1391 ) has a singularity at p = ^Ilb. This is the 
location of the horizon at the zeroth order and the zeroth order metric itself is not 
invertible here. 

The first order perturbation has a log piece which also blows up here. This 
singularity could be just a coordinate singularity in which case it could be re- 
moved by going to a different coordinate system as it happened for the boosted 
black brane, or it could be a real singularity. If it is a real singularity, it is naked 
because it coincides with the original horizon at late time. At late times the solu- 
tion approaches a boosted black brane but since the horizon coincides with a real 
singularity, no infalling observer can continue life after reaching the horizon. 

To analyse the singularity in the Fefferman-Graham coordinates we will sim- 
ply translate the metric to Eddington-Finkelstein coordinates (r,x). It will be of 
course suffice to change our coordinates near p = V2Z?, however, for the sake 
of completeness and better general understanding we will do the change of co- 
ordinates exactly upto first order in the derivative expansion. The Eddington- 
Finkelstein metric which we will get as a result of this translation will also be 
an exact solution of Einstein's equation upto first order in x-derivatives. We now 
return to the equations (fT5l) in section 3 which gives the translation between the 
two coordinate systems. We still treat the Fefferman-Graham coordinates (p(x, r), 
Z v (x, r)) as unknowns, but the third unknown is now the G MV (x, r) which appears in 
the Eddington-Finkelstein metric (fTTT) . The zeroth order solutions to these three 
are known and are given in ([T9l and (1241) . To find the corrected solutions due to 
change in the Fefferman-Graham metric at first order it is straightforward to per- 
turb these equations and solve them exactly at first order. The complete solutions 
to the three unknowns exact upto first order are: 



P = 
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t = x? + u"bk(br) + u^b 2 k A (br) + (u.d)u"b 2 k B (br) (44) 



1 2 

G^ v =r 2 P MV + {-r 2 + -g^u^Uy + l/bFibf)^^ - r((u.d)(u M u v ) - -u^Uyid.u)) 



(7-1)6 2 1 

where, 



1 x + 1 

= -{log{ ) - 2arctan(x) + n) (45) 

4 x - 1 

, Kj ,(x+i)V + i) 

r (jc) = -{log{ ) - Larctanyx) + n) 

4 x 4 

and £a(x), feW satisfy the following differential equations 

^ = --^-AKx) + -f=) (46) 
ax xr - 1 Yx 4 - 1 

Jfce 1 k{x)x 2 

dx xVx 4 - 1 - 1 

with the boundary condition that they vanish at x = oo. One may easily check that 
if we do the Taylor series expansion of p,z M in 1/r, we can reproduce the results 
(|29l) of section 3 in which we have solved these equations using a power series 
ansatz. 

The crucial point, as realized by authors of [10] is that in the Eddington- 
Finkelstein coordinates if there is a blow-up in G MV (x, r) it should be a real sin- 
gularity. For a general conformal fluid at first order with rj/s = y/An, the cor- 
responding solution in Eddington-Finkelstein coordinates has G /jv (jc, r) given by 
(14"4"1) . Except for the log term which appears in the last line, all other terms are 
well behaved for r > and the log term blows up at r = 1 /b, the location of the 
unperturbed black brane horizon. Only when 7=1, the coefficient of the log term 
vanishes and so the naked singularity at r = l/b is absent. For this value of y we 
have in fact reproduced the of the Eddington-Finkelstein metric given by the 
authors of [10]. 

We learn the following general facts. The translation to Eddington-Finkelstein 
coordinates exists for an arbitrary solution in the Fefferman- Graham coordinates 
irrespective of whether there is any naked singularity or not. Also the Fefferman- 
Graham coordinates have a power series expansion in terms of the inverse of the 
radial Eddington-Finkelstein coordinates for all cases. For all cases, the change 
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of coordinates also become singular at the location of the original horizon in the 
Eddington-Finkelstein coordinates which is r = 1 /b. 

We can continue the regularity analysis to higher orders in the derivative 
expansion by solving the equations (TT5T) for translating the solution from the 
Fefferman-Graham coordinates to Eddington-Finkelstein coordinates order by or- 
der in the derivative expansion as well. In this way at each order we will be able 
to determine what values the coefficients in the terms of the hydrodynamic stress 
tensor should have so that a naked singularity is avoided. It would be interesting 
to see if we can understand the values of these coefficients of the hydrodynamic 
stress tensor, more directly in terms of the geometry of the unperturbed boosted 
black brane horizon. 

We will conclude this section by emphasizing certain points. 

• We can think of translating to outgoing Eddington-Finkelstein coordinates 
also as an attempt to remove the singularity and then as expected the situa- 
tion will be time-reversed. We will now need y = -1 for regularity. In the 
boundary theory, all fluid dynamical solutions will then be time-reversed 
and our gravity solutions will be perturbed white-hole solutions exact upto 
first order in the derivative expansion. 

• We could have attempted to fix y by studying regularity at the horizon by 
computing curvature invariants (like R Mvpo -R Mypcr ). However, we do not know, 
if for these "hydrodynamic" space-times, checking that a finite number of 
curvature invariants do not blow up at the horizon will suffice to demon- 
strate regularity. So the best strategy is to translate to a coordinate system 
where the solution is explicitly regular upto first order in the derivative ex- 
pansion and this is what we have done here. For the sake of completeness, 
however, we have studied a few curvature invariants and have found that the 
leading singularity of R tlVfxr R , * vp(r at second derivative order vanishes for the 
right choices of y which are 1 and -1, the details of which are presented in 
Appendix B. 

• The derivative expansion in Fefferman-Graham coordinates is equivalent to 
the same in Eddington-Finkelstein coordinates to all orders in the derivative 
expansion even when the solutions do not have a regular horizon. This is so 
because the equations (fT5l) for translating Fefferman-Graham coordinates to 
Eddington-Finkelstein coordinates can always be solved order by order in 
the derivative expansion as well. In fact, this is natural, because any asymp- 
totic AdS solution can be written in the Fefferman Graham coordinates. 
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6 Discussion 



We will point out some implications of our results for dynamics in the universal 
sector of CFT. Our first result is that a solution of pure classical gravity is uniquely 
specified by the stress tensor. This implies that the dynamics of all states in the 
universal sector of the dual CFT at strong coupling and large N is completely 
determined by the conservation of the traceless stress tensor. The implication for 
dynamics on the CFT side is even more surprising than the result for classical 
theory of gravity itself. It is surprising because to characterise a state uniquely 
we would typically need the expectation values of infinite number of operators. 
However, it is not hard to give an example of a special sector of states with this 
property in a 2D CFT. These special states are spanned by L n \VAC > in ^ 2)and 
are created by descendendants of the identity operator (L„),with n > 2, acting on 
the vacuum. Each such state is uniquely character by the L eigenvalue n, hence 
by the expectation value of the stress tensor T(z). Moreover each state LJVAC > 
(n > 2) being an eigenvector of the Hamiltonian, the sector spanned by these states 
is closed under time evolution. It would be interesting to find such examples 
of class of states in CFTs in higher dimensions where the expectation value of 
the stress tensor uniquely identifies each member and moreover is closed under 
time evolution. The real question, however is, whether we can give an intrinsic 
microscopic description of the universal sector of CFTs with gravity duals. If we 
can achieve this, we will be able to understand better how the vev of the stress 
tensor and its conservation alone determines the dynamics in the universal sector 
completely. 

Our second set of results are (a) all hydrodynamic stress tensors are free of 
the "abed" type of pathology, which means that they preserve the asymptotic AdS 
boundary condition and (b) there is a unique hydrodynamic stress tensor for which 
there is no naked singularity. This means that the late time equilibriation in the 
boundary CFT can be determined by a unique and universal hydrodynamic stress 
tensor. The coefficients of the terms should be set to values which avoids for- 
mation of naked singularities in the bulk. It would be interesting to find out an 
intrinsic microscopic definition for the higher order coefficients of the hydrody- 
namic stress tensor, in terms of say, multi point correlations of the stress tensor. 
The first order coefficient, namely the viscosity has indeed such a definition in 
terms of two-point correlation function of the stress tensor and the validity of the 
definition can be verified by the AdS/CFT correspondence as well. So we may 
hope that a pure gravity analysis should suffice to arrive at similar definitions for 
the higher order coefficients in the hydrodynamic stress tensor. 
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We would like to mention that while we were updating our work, it was a 
great pleasure to find out that our method has been generalised in [20] to com- 
pute the stress tensor in the universal hydrodynamic sector of strongly coupled 
large N dual theories of various p-branes, which are in most cases non-conformal. 
We would like to take this opportunity to mention that since our method keeps 
the asymptotic boundary condition manifest, it could be given a preference when- 
ever implementing the asymptotic boundary condition in Eddington-Finkelstein 
coordinates becomes laborious or complicated. 

Finally, we would like to point out, that it will be interesting to find a physical 
understanding of the "abed" type of pathology. Our results in the hydrodynamic 
sector gives support for claiming that whenever we have late-time equilibriation 
in the boundary stress tensor, this pathogy is absent. It will be interesting to find 
a real example with such a pathology and trace its physical origin. 
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Appendix A : Proof of the power series solution for 
AdSs asymptotics 

Here we will prove that any asymptotically AdS 5 solution of Einstein's equation 
with a negative cosmological constant, in the Fefferman- Graham coordinates, has 
a solution for g^ v which is a power series in the radial coordinate when the bound- 
ary metric is flat. Though not explicitly mentioned in most of what follows, it 
should be kept in mind that here we are specifically investigating five-dimensional 
solutions with a flat boundary metric. At the end, we will mention if our proof can 
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be generalised to other cases. 

To simplify the proof we first rearrange the tensor and the scalar components 
of Einstein's equation © while keeping the vector components of Einstein's equa- 
tion unchanged. The old scalar equation is added with an appropriate linear com- 
bination of the trace of the old tensor equation so that now it does not contain any 
term which has second derivative of g^ v with respect to the radial coordinate p. 
Since the vector equation also does not contain any term with second derivative 
of g M y with respect to the radial coordinate we can now think of the vector and 
scalar components as a set of five constraint equations. We also change the tensor 
components of Einstein's equation by appropriately replacing Tr(g~ l g') using the 
new scalar equation. We do this so that now the tensor equation by itself is suf- 
ficient to determine all the p" coefficients of g^. The old tensor equation had the 
feature that to determine g^ v , the coefficient of p 8 in g^, we had to use the scalar 
equation as well, but now this can be fully determined using the tensor equation 
alone. So our equations now are as below. 



It is not difficult to see that we can use a power series ansatz to solve the tensor 
equation as at the n-th order. At the n-th order the only terms which can contain 
g ( „ v , or Tr{g (n) )T]^ are gfa,, g' MV and Trig^g^g^. Now since the tensor equation 
contains no term with Tr(g~ l g')g,j V , at the n-th order,for n > 4, the tensor equation 
gives us n(n - 4)g (n)M J2 = f {n)tlv {t po ), where f (n)MV (t pcr ) is a polynomial in t p(T and 
its various derivatives with respect to the boundary coordinates only. Hence, for 
n > 4, we can always solve g( n ) MV using the tensor equation alone. 

We have now got to show that the power series we have so obtained as a 
solution to the tensor equation is consistent with the vector and scalar constraints. 
We will do this by the method of induction iterating over the various coefficients 
of p n in g^y, order by order in n. We will first establish the following fact that the 
p-derivative of the vector and scalar constraints vanish when the tensor equation 
along with the vector and scalar constraints are satisfied. This just articulates 
the intuition that once the initial data consisting of g pv and g' satisfy the vector 



\g" ~ ^~g' ~ ^ V + jTr(g- l g')g' - Ric(g) 
2 2p 2 4 

+ gi^R(g) + ^Tr{g- l g'g' l g ') - _L(7>ory )) 2 ] = o 



(47) 



V>ory)-vx, v = o 

R{g) + -Tr(g- l g') + jTr(g~ l g'g' l g') - l[Tr(g' l g')] 2 = 
p 4 4 



(48) 



(49) 
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and scalar constraints on hypersurface with a fixed value of the radial coordinate 
p, the dynamical evolution in p should be such that the constraints should be 
automatically satisfied for any other hypersurface. To show this we will need the 
following: 

Kr' = jgTFvgL + v^ v - v^) (50) 

One can use the tensor (1471) and scalar (|49l) equations to write 

— ~Z&yR = Sav ~ 7T~ 8^ §av ~ 7^8** SafiS^Syv + T-^ r C? 8 8av 

9 (51) 

+ 4-Tr(g- l g'W v - l -K\Tr{g- l g") - Tr(g~ l g' g~ l g') + \{T r{g~ l g')f] 

Now when all the equations (1471) . ([481) and (1491 ) are satisfied, the p-derivative of 
the vector constraint can also be written as: 

(V^TVOf V) - V^v)' =d,[Tr(g- i g" - jg"V«"V) + \{Tr(g- V)) 2 ] 

4 7 (52) 

" Vv(g av g'; a - g ap g'p y g" v g' a , + 2S Va 8' m Tr(g- l g')) 

Now comparing the right hand sides of (I5TT) and (|52l ) using all the equations of 
motion again, we see that 

(V,Tr(g- l g') - V v g'J = V V (R; - \s;R) (53) 

So the Bianchi identity implies that the p-derivative of the vector equation should 
vanish when all the equations of motion are satisfied. We will now get to the scalar 
equation. 

When the vector equation of motion (|48l) is satisfied we get 

V = -\R m (g aP g'p v ) + \KJ£*fa) + \v^ v Tr{g- l g') - Vs; v (54) 
This implies that when the vector equation of motion is satisfied, we have: 

R' = -g" v g' v X a R, a (55) 
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On the other hand the vanishing of the p-derivative of the scalar constraint ( |49|) 
ought to give us: 

R'=~ \Tr{g- l g'g- l g") + l-Tr{g- l g'g~ l g') 

+ \Tr{g~ l g'g~ l g'g- l g') - l -Tr{g~'g')Tr{g- l g'g-'g') + ^-[Trig' 1 g')f 
2 4 2p 

(56) 

Now using the tensor and scalar equations of motion, we can see that the right 
hand sides of (1551 ) and (|56l ) are the same, or in other words the p-derivative of the 
scalar constraint indeed vanishes when all the equations of motion are satisfied. 
So we have established that the p-derivatives of the all the five constraints vanish 
when all the equations of motion are satisfied, or to state compactly 

(ST]), (H8), dH => W, W (57) 

To prove that the power series solution of the tensor equation is consistent with 
the constraints, we will use the above at p = 0. To obtain a condition for g^ v 
(the coefficient of p" in g MV ) from the tensor equation we need to differentiate it 
n-2 times with respect to p and then set p = 0. Similarly to obtain a condition for 
g(n)fiv from the vector and scalar constraints we need to differentiate each of them 
n-1 times with respect to p and then set p = 0. 

The vector and scalar constraints imply that g(2 )M v should vanish while the ten- 
sor equation identically vanishes at this order. The tensor equation for g(4) flv which 
we have appropriately renamed t MV , also identically vanishes while the vector con- 
straint gives us the conservation equation <9% v = and the scalar constraint gives 
the tracelessness condition Tr(t) = 0. We can start our induction from here, since 
the three equations are all consistent with each other upto this order 

Let us suppose, by the induction hypothesis that the solution for g(„-\ )llv ob- 
tained from the tensor equation is consistent with the vector and scalar constraints. 
We now denote the m-th p-derivative as m . So, by induction hypothesis, the three 
equations (n - 3Y(|47Kp = 0), (n - 2)'<HEKp = 0) and (n - 2)'<g§Kp = 0) are 
consistent with each other. Now we iterate by determining g^v from the tensor 
equation, or in other words we solve 

in - 2)'gj}(p = 0) (58) 

But by induction hypothesis we can assume (n-2) / (|49l)(p = 0) and (n-2) / (|48l(p = 
0) are consistent with the tensor equation. Now our result (1571 ) for a general fixed 
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p hypersurface implies that 

(n - 2)'<@7J, (n - 2)'<m, (n - 2)'@9) => (n - l)'®, (n - 1)'© (59) 

We can apply the above at p = 00 to iterate and say that if the solution for g( n -i)nv 
from the tensor equation is consistent with the constraints so would the solution 
for g(n)nv from the tensor equation be. This completes the proof by induction that 
the power series solution of the tensor equation is consistent with the constraints. 

Let us see if our proof can be generalised to other cases, in particularly for all 
dimensions if the boundary metric is flat. The only change in the equation of mo- 
tion happens to be the coeffiecient of g' MV in the tensor equation. Let us, for exam- 
ple, take the case when the number of boundary coordinates is six. We can check 
by hand that all g^v vanish for all n such that < n < 6 and g^v cannot be deter- 
mined from the tensor equation for an exactly similar reason as for g^y when the 
number of boundary coordinates was four, namely the tensor equation identically 
vanishes. The vector and scalar constraints imply conservation and tracelessness 
°f g(6)/jv implying that it should be identified with the stress tensor (and indeed 
it has been shown in [5] that this agrees with with the Balasubramanian-Krauss 
stress tensor). We can begin our induction, from here as before and hence our 
proof generalises. So, the general problem in applying the induction is to show 
that the equations of motion are consistent with the power series ansatz at g^v 
We have not been able to prove it generally but we have checked it upto d = 6. 
The same problem appears when we try to apply induction to prove the validity 
of the power series solution when the number of boundary coordinates is odd, but 
the boundary metric is arbitrary. Before we apply induction, we need to prove 
that the power series works at g(d) flv , (in fact this is harder to show, because when 
the boundary metric is not flat g^yS do not vanish for < n < d). However, 
Fefferman and Graham have proved the validity of the power series solution by a 
different method for an arbitrary boundary metric when the number of boundary 
coordinates is odd. 

9 At p — the statement (T59b has a non-trivial content strictly for n > 2, because of the 
slight technicality that what we really need to use to find a condition for g( n )^v is that we need 
to differentiate (jO(|49l>) not really d49l l n-1 times. So at p — 0, this result is trivial for the scalar 
constraint when n — 2 and we do not need to use the result d59l l, but since the first step of induction 
starts from n — 4, it is safe to use this in the iteration procedure. 
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Appendix B: On fixing rj/s by calculating curvature 
invariants 



We have already done the regularity analysis of our first order solution in Fef- 
ferman Graham coordinates by translating to Eddington-Finkelstein coordinates 
where the regularity or irregularity becomes manifest. However, one may ask if 
the regularity analysis can be done also by calculating some curvature invariants. 
We will see that indeed at the first order, this analysis can also be done by cal- 
culating an appropriate curvature invariant, but we will argue that there may not 
be a finite number of curvature invariants which can be reliably used to fix all the 
coefficients in the hydrodynamic stress tensor at higher orders in the derivative 
expansion. 

Before we do that, we want to point out that though the metric in Fefferman- 
Graham coordinates and in Eddington-Finkelstein coordinates could be made co- 
ordinate equivalent upto any given order in the derivative expansion for an arbi- 
trary hydrodynamic stress tensor, the curvature invariants calculated from the two 
metrics will typically never be the same! Let us examine why this should happen 
at the first order itself. Any typical curvature invariant, like the Ricci scalar R it- 
self, will show a divergence only when we expand it to second order in derivatives 
of the boundary coordinates. In this case, this should be so, because the metric in 
either coordinate system is a solution of the equations of motion upto first order 
in derivatives of boundary coordinates. However, the second order piece in R cal- 
culated from the metric in either coordinate system will not be the same, because 
the two metrics are related by a coordinate transformation only upto first order 
in derivatives. In fact we will explicitly demonstrate that R itself can be used to 
fix the value of 4ttt]/s in the Eddington-Finkelstein metric at first order but not 
in the Fefferman- Graham metric at first order. So the procedures of using curva- 
ture invariants to fix the coefficients in the hydrodynamic stress tensor in the two 
coordinate systems are indeed very different! 

Another crucial aspect should be kept in mind because this also features in 
comparing curvature invariants calculated from the metrics in the two coordinate 
systems. Fundamentally, solving Einstein's equations in either of the two co- 
ordinate systems involves a trade-off between manifest regularity and manifest 
asymptotic boundary condition. The solution in Eddington-Finkelstein coordi- 
nate system at the zeroth order and also at the first order for the right value of 
4nr]/s are manifestly regular so any curvature invariant calculated at the hori- 
zon will be regular to all orders as well. However, the solution preserves the 



34 



asymptotic AdS boundary condition only upto first order in derivatives as it can 
be translated to Fefferman-Graham coordinate system only upto that order. The 
solution in Fefferman-Graham coordinate system at first order, of course preserves 
boundary condition to all orders, but even for the right choices of 4nr]/s it is not 
regular to all orders. In other words, for the right choice of 4nrj/s all order diver- 
gences should vanish when we calculate curvature invariants from the metric in 
Eddington-Finkelstein coordinate system, but in case of the solution in Fefferman- 
Graham coordinates at first order, at most the leading divergence at the second 
order vanishes for the right choice of 4nr]/s. In fact, for certain curvature invari- 
ants even that do not happen. Of course, eventually if we add a right second order 
correction to the Fefferman-Graham metric, all divergences in the curvature in- 
variants at the second order should vanish, but still divergences at higher orders 
will remain and so on. We will illustrate the first order case with examples below. 

To compute curvature invariants it is useful to first choose a velocity and tem- 
perature profile. As mentioned before, the vector constraint in Einstein's equations 
of motion demand that the velocity-temperature profile should be a solution of the 
relativistic Euler equation 

d M b d.u 

— = (u.d)u M -u M — (60) 
b J 

We call our boundary coordinates (t, x, y, z) and we select the following static 
velocity profile which is a relativistic version of laminar flow 

if= . 1 (l,gy,0,0) (61) 
VI - a 2 y 2 

where a is a constant of dimension 1/length. The advantages of using this velocity 
profile are twofold, namely, 

• The relativistic Euler equation gives us that temperature, hence b, should be 
a constant. 

• It is easy to employ the derivative expansion by using the following trick. 
We note that the only non-trivial derivatives of the boudary coordinates are 
the y-derivatives. Any y-derivative of the velocities will bring in an extra a 
which is unpaired with a y so that it picks up the right dimension. Hence to 
do the derivative expansion we may first set y = p/a and simply do a Taylor 
expansion in a about a = 0. The correct dimensionless parameter of the 
derivative expansion, of course will be ab. 
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We can use the above velocity-temperature profile in the first order solution in any 
coordinate system. Though away from the boundary the boundary coordinates 
(or, in other words, the field-theoretic coordinates) in a given coordinate system 
will mix with all the coordinates in another coordinate system, at the boundary 
they will always align with other. This is, how solutions in two different coor- 
dinate systems come to share the same boundary stress tensor and also the same 
conservation equation, which in this case, is the relativistic Euler equation. 

If we use the above velocity-temperature profile to calculate R in the Edding- 
ton Finkelstein coordinate system we will find that 

R = -20 + a 2 



8(1 - a 2 y 2 ) 2 b 4 r 6 

(y - l)(fr V(9 + 3y - In) - \6b 5 r 5 + 2nb 6 r 6 )) + ^ 
(br - 1)(1 + br + b 2 r 2 + b 3 r 3 ) 

(y - l)(y + 1 - Sb 3 r 3 )b 2 r 2 Log(l - + 0(1)] + 0(a 3 ) 

At the zeroth order in a, R should of couser be -20 and at order a, R should of 
course vanish because our metric is a solution of equations of motion upto first 
order. At order a 2 , we indeed expect some divergence at the horizon, which is at 
r = I /b, because the metric is explicitly not regular there unless y = Anrjl s = 1. 
We see that when 4m]/ s = y = 1 all divergences go away. This feature replicates 
also at higher orders in a. On the other hand, if we calculate R from the 
Fefferman-Graham metric at first order, we get 

2 \2%b w p\\2b 4 y 2 + 4b 2 p 2 + 3yV) 

+ Cl[ (1- fl 2 y 2 )W - p4)2 (4 £4 + p4)3 

16&W 4b 4 -p 4 , 1 J 

+ (l-aW(4M + p^ L ° g( 4MT7 )] + ° (a } 

At order a 2 , we see that there is a leading inverse power two divergence for any 
value of y and a subleading log divergence except when y = 0. So this is useless 
to figure out the right value of y. Of course this will certainly be useful to fix 
certain coefficients of the hydrodynamic stress tensor at second order, because 
these divergences should go away for any right second order correction to the 
Fefferman-Graham metric. 



10 We would like to thank Sayantani Bhattacharya for confirming that this indeed happens for 
arbitrary velocity and temperature profiles. 
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It turns out, however, that, R fiyp(r R IJVp(r can be used to fix the value of y in the 
Fefferman- Graham metric. We get 

^ w _ 4(1280ft 16 + m0b i2 p 4 + 2784fr 8 p 8 + 80ft 4 p 12 + 5p 16 ) 
R ^ R ~ (4M+p 4 ) 4 

, 2{\-y 2 )b 6 1 /- 
- a\ ) Y) + 0( — ) + 0(Log( y(2b - p)) + 0(1)] 

(1 - a 2 y 2 ) 2 (p - y[2b) A (p - ^W) 2 

(64) 

We see that the zeroth order piece is always finite and independent of y and at 
order a (for some reason we do not understand) the scalar vanishes. However, 
at order a 2 , we find that when y is 1 or -1 the leading divergence at p = ^llb 
goes away, though, the subleading divergences remain and as before, they should 
disappear when we add any right second order contribution to the Fefferman- 
Graham metric. We are also not sure, if by computing R Mvp(T R Myp ' T itself we can fix 
the values of all the coefficients in the hydrodynamic stress tensor at second order. 
To fix all the coefficients of the second order hydrodynamic stress tensor, one may 
have to look for another appropriate curvature invariant. 

It is certainly, worth exploring, if the "hydrodynamic" Fefferman- Graham so- 
lutions are '"special" enough so that computing a finite number of curvature in- 
variants will suffice to determine regularity, hence in fixing all the coefficients 
in the hydrodynamic stress tensor to all orders. We will leave this for a future 
work. Nevertheless, our procedure of fixing the coefficients in the hydrodynamic 
stress tensor by translating to Eddington-Finkelstein coordinate system works for 
all orders in the derivative expansion. 
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